(Meta-) stable reconstructions of the diamond (111) surface: interplay between 

diamond- and graphite-like bonding 
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Off-lattice Grand Canonical Monte Carlo simulations of the clean diamond (111) surface, based 
on the effective many-body Brenner potential, yield the (2 x I) Pandey reconstruction in agreement 
with ab-initio calculations and predict the existence of new meta-stable states, very near in energy, 
with all surface atoms in three-fold graphite-like bonding. We believe that the long-standing debate 
on the structural and electronic properties of this surface could be solved by considering this type 
of carbon-specific configurations. 



C3 



-a 
o 

> 
On 

m 

O 
O 
O 

-1— > 



1 

G 

O 

o 



X 



The discovery of fullcrcnc has awakened an increasing 
interest in carbon based nanostructures,as well as in pro- 
cesses, such as diamond graphitizatiorfl or the graphite- 
to-diamond transformationB leading to structures, which 
promise to have desirable properties of both graphite and 
diamond. It is important to develop predictive schemes 
to treat diamond, graphite and mixed bonding with ap- 
proaches able to deal with large structures, often beyond 
the possibility of ab-initio calculations. The effective 
many-body empirical potentials due to Tersofm for group 
IV elements are very accurate for Si and Ge, but less reli- 
able for C. In particular, for C the Tersoff potential yields 
the unreconstructed (111) (1 x 1) surfaca-ias the most 
stable against the expejrimental evidenceoEl of a (2 x 1) 
Pandey reconstruction!]'!! For the (001) face it predicts 
dimcrization with a strong asymmetric displacement of 
the third-layer atoms .El Here, we use the potential pro- 
posed by BrennerEJ (parametxi^ation I) and show that it 
is reliable also at the surface £3 

We perform off-lattice Grand Canonical Monte Carlo 
(GCMC) simulations!^ of the (111) surface of diamond. 
We find the unbuckled undimerized (2x1) Pandey chain 
reconstruction as the minimum energy structure and 
three new mcta-stable states, close in energy, with all sur- 
face atoms in a three-fold graphite-like bonding. Two of 
them are obtained by a strong dimcrization of the lower 
(4-fold coordinated) chain, inducing a small dimerization 
of the upper (7r-bonded) chain as well. The third meta- 
stable (V3 x V3)R30 o reconstruction is formed by a reg- 
ular array of vacancies. 

Surprisingly, the reconstruction of clean diamond(lll) 
is not yet established in detail, although there is a consen-i 
sus that the 7r-bonded Pandey (2 x 1) reconstructiontm 
is the most stable. One important issue is whether 
this surface— is—metallic or semiconducting. In most 
calculationsLrEZl the band of surface states is metal- 
lic whereas experimentally the highest i©«pipied state 
is at least 0.5 eV below the Fermi level.t§ll3 Dimeriza- 
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ang the 7r-bonded chain could open the surface 
gapE^tcil but only one total-energy calculation— obtains 
slightly dimerized chains yielding a 0.3 eV gapE^I in the 
surface band. Experimentally, recent X-ray dataQ and 
medium-energy ion scattering^ do not show any dimeriza- 
tion but favor the (2x1) reconstruction accompanied by 
a strong tilt of the 7r-bonded chains, similar to the (2x1) 
reconstruction of Si(lll) and Ge(lll). The,tiltjs how- 
ever not confirmed by theoretical studies .0~tZft3 Also, 
relaxations in deeper layers are debated. The bonds be- 
tween first and second bilayers are foundjio be elor 
by an amount which varies between l9cEn3 and S 
Bonds between the second and third bilayers are slightly 
shifted (< 1%) in theoretical studies while X-ray data 
suggest a 5-6% relaxation.Qtl 

Experimentally, uncertainties can be caused by vari- 
ations—of. surface preparation. A pastial graphitiza- 
tioiJJl'tjO and other structural phasesEiaEa can coexist 
at the real surface. It is noteworthy that .m o s t structural 
models were first suggested for Si and Ga3>E3Ea and then 
extended to diamond.QO However, the former always fa- 
vor tetrahedral four-fold coordination whereas C favors 
also the graphite-like three-fold bonding, the latter being 
in fact energetically stable in the bulk at normal condi- 
tions. One can therefore expect diamond to have addi- 
tional low-energy surface structures, such as the meta- 
stable reconstructions presented here. 

Empirical potentials, although less accurate than ab- 
initio calculations, allow to explore larger portions of 
phase space and can lead to unexpected structures, which 
can then be tested in more accurate calculations and 
taken into account in the experimental data analysis. We 
exploit the MC simulated annealing scheme to overcome 
potential energy barriers and identify low-energy surface 
structures. Atoms in the bottom layers are kept fixed 
at their ideal bulk positions while the others (usually, 
four bilayers of 128 atoms each) are mobile. We con- 
sider either a VNT, PNT or PfiT ensembles for differ- 
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ent tasks. The canonical VNT ensemble is used to mini- 
mize ordered structures. During the simulated annealing 
we allow for volume fluctuations (PACT). We consider 
also a grand canonical P/iT ensembleEJ to access struc- 
tures with different number of atoms than in the bulk 
terminated structure. Each atom creation/destruction 
is enabled only in the near-surface region (2-3 top bi- 
layers) and followed by 1000 MC equilibration moves. 
We have optimized the implementation of the potential 
by use of neighbor lists which allow to calculate energy 
variations on a finite portion of the sample. The one- 
dimensional functions defining the potential are stored 
in tables with a fine grid and calculated by linear in- 
terpolation. The attractive Va, repulsive Vr and the 
cut-off f c terms are stored as a function of the square 
of interatomic distances to avoid square root operations. 
The three-dimensional F function is. stored on a finer 
grid and the tricubic interpolation!^ is replaced by a 
linear one, reducing the terms to be evaluated from 64 
to 8. We note that, in a strongly covalent system like 
diamond, creation and destruction are very improbable 
events, also because immediately after creation or de- 
struction the neighboring atoms have not yet adjusted to 
the new local environment; after a destruction the sys- 
tem can gain up to 2 eV by-jrelaxation. This energy has 
been added as an umbrellaEH in the acceptance rule for 
creation/destruction. 
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FIG. 1. Left: relaxed lxl, right: Pandey 7r-bonded chain. 
The top (111) and side (110) views are presented in the top 
and bottom images, respectively. The relative changes of the 
bond lengths and the atom labeling schemes used in Table [j] 
are indicated. 

Throughout, we give energy gains AE per lxl unit 
cell, relative to the bulk-terminated surface. We find 
the relaxed (lxl) and Pandey (2x1) structures shown 
in Fig. H to have AE = 0.244 and 1.102 eV, respec- 
tively. Apart from a 4% elongatiorLaLthje bond between 
first and second bilayer against 89cDOa for the Pandey 
structure, o ur resuL ts-agree remarkably well with ab-initio 
calculationsMliaMEl 

To the best of our knowledge, no one has succeeded 
before in simulating a spontaneous transition from the 



ideal (1 x 1) to the (2 x 1) reconstructed diamond (111) 
surface. The (2 x 1) reconstruction of both diamond(lll) 
and Si(lll) is associated with a large coherent displace- 
ment of the first bilayer by more than 0.5 A, accompanied 
by rebonding in which one atom in the 2x1 surface unit 
cell changes coordination from four to three while an- 
other does_the opposite. Ab-initio molecular dynamics 
simulation^] of the spontaneous (2x1) reconstruction of 
Si(lll) shows that breaking and formation of a new bond 
occurs at the same time. We find that for diamond, in- 
stead, the bond breaking leading to an increase of three-, 
fold (graphite-like) atoms precedes the bond formation.E2l 
In such a situation, competition between the (2 x 1) re- 
construction and-ipartial) surface graphitization becomes 
very important ,E3 especially if the annealing is performed 
at high temperatures. Conversely, a low annealing tem- 
perature makes it very difficult to overcome the potential 
barrier between the two ordered structures. In Fig. ^| we 
show the top view of a diamond (111) sample, obtained 
from the ideal relaxed (lxl) structure after an anneal- 
ing cycle (about 0.5 • 10 6 MC steps) at T = 750 K. The 
efficiency of phase,space exploring is improved by increas- 
ing the step sizeO as to reduce the acceptance rate from 
the usual 50% down to 20-25%. In Fig. || we emphasize 
pieces of the lower, four-fold coordinated Pandey chains 
which represent the final stage of the (1 X 1) — > (2 X 1) 
transition. Only one rotational domain is present and 
the relative position of the formed chains is somewhat 
disordered. 




FIG. 2. Top view of a reconstructed diamond(lll) sam- 
ple after a spontaneous transition from the bulk-terminated 
(lxl) structure. The four- fold coordinated atoms in the first 
bilayer (atoms 13 and 14 in Fig. [l]) are shown as dark balls and 
connected by dark solid lines. Periodic boundary conditions 
within the surface plane are used. 
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FIG. 3. Meta-stable structures of the diamond(lll) surface arranged as in Fig. Left and central structures are the (2 x 1) 
and (4 x 1) dimerized reconstructions. Short and long distances in the dimerized chain are indicated by solid and broken lines, 
respectively. The structure on the right displays the vacancy (V3 x VS)R30° reconstruction. 



TABLE I. Ceerdinates of atoms (in A) within top three bilayers for the Pandey model, (2 x 1) and (4 x 1) dimerized 
reconstructionsESI The Brenner potential gives the diamond bulk bond length of 1.541 A and the 2x1 unit cell has in-plane 
dimensions 4.357 x 2.517 A 2 . The numbers in brackets for the (4 x 1) dimerized reconstruction are coordinates of atoms in the 
second half of the unit cell (denoted with primes in Fig. []). 
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11 


3.064 


-0.943 


5.525 


3.071 


-0.942 


5.638 


3.072 (7.435) 


-0.532 (-1.355) 


5.648 


(5.645) 


12 


2.370 


0.315 


5.519 


2.353 


0.332 


5.625 


2.379 (6.742) 


0.728 (-0.095) 


5.646 


(5.641) 


13 


1.007 


0.313 


4.792 


1.177 


0.594 


4.804 


1.210 (5.591) 


0.714 (-0.081) 


4.802 


(4.793) 


14 


0.076 


-0.941 


4.804 


-0.112 


-1.188 


4.820 


-0.151 (4.228) 


-1.325 (-0.561) 


4.813 


(4.806) 


21 


3.816 


-0.942 


3.325 


3.835 


-0.965 


3.331 


3.838 (8.202) 


-0.911 (-0.989) 


3.330 


(3.325) 


22 


3.041 


0.316 


2.880 


3.065 


0.299 


2.872 


3.072 (7.436) 


0.350 ( 0.275) 


2.866 


(2.862) 


23 


1.557 


0.314 


3.290 


1.573 


0.341 


3.297 


1.589 (5.967) 


0.378 ( 0.250) 


3.313 


(3.302) 


24 


0.895 


-0.943 


2.737 


0.896 


-0.910 


2.759 


0.898 (5.272) 


-0.878 (-1.015) 


2.774 


(2.767) 


31 


3.065 


0.316 


1.300 


3.075 


0.318 


1.295 


3.084 (7.449) 


0.322 ( 0.306) 


1.299 


(1.291) 


32 


2.346 


-0.942 


0.760 


2.251 


-0.946 


0.770 


2.361 (6.727) 


-0.938 (-0.955) 


0.769 


(0.765) 


33 


0.883 


-0.944 


1.227 


0.889 


-0.930 


1.240 


0.895 (5.265) 


-0.918 (-0.974) 


1.252 


(1.246) 


34 


0.152 


0.315 


0.739 


0.156 


0.328 


0.752 


0.162 (4.533) 


0.342 ( 0.293) 


0.756 


(0.756) 



At higher temperatures, apart from a tendency to- 
wards graphite- like structures, we observe also other 
ordered phases, similar to the Pandey reconstruction 
but accompanied by a strong dimerization of the lower 
atomic chainEj This dimerization can be performed in 
two ways, leading to the (2 x 1) and (4 x 1) structures 
shown in Fig. |^. The atom coordinates are givenE3 in Ta- 
ble | along with those for the Pandey reconstruction. The 
energy gain AE is 0.883 and 1.023 eV, for the metastable 
(2x1) and (4x1) respectively. In the dimerized chains 
the short/long distances between atoms 13 and 14 are 
1.459/2.215 and 1.444/2.455 A for the (2 x 1) and (4 x 1) 
instead of 1.562 A in the Pandey reconstruction. This 
dimerization induces dimerization of the 7r-bonded chain 
as well, albeit small (< 1%), which might be of impor- 
tance for surface electronic properties. The (4x 1) dimer- 
ized structure is only 160 meV (« 2000 K) per broken 
bond higher than the Pandey structure, i.e. these phases 
can coexist at the surface at high temperatures. Indeed, 
by heating the ordered Pandey structure to 2350 K we 
observe a partial transformation to the dimerized state as 



shown in Fig. [4|. Note also a precursor of surface graphi- 
tization. 

Lastly, in GCMC runs we find the (V3 x ^#30° 
reconstruction formed by an ordered array of vacancies 
as shown on the right hand side in Fig. |[ The bond 
lengths between the atoms in the first bilayer is reduced 
to 1.390 A, slightly less than the equilibrium bond length 
in graphite (1.42 A). The bonds between the first and sec- 
ond layer are elongated by ~ 1% while the other bond 
lengths are close to the bulk value. Taking the bulk bind- 
ing energy 7.346 eV as the chemical potential, the energy 
gain is estimated to be AE = 0.6145 eV. Once formed, 
this structure is found to be very stable and remains un- 
altered after long GCMC annealing cycles at T = 
K.Ell Similar structures have been discussed for SiEiro 
For diamond(lll) a (2 x 2) vacancy structure was shown 
to be energeticallyr-unfavorable compared to the relaxed 
(1 x 1) structure.^ However, contrary to the (2 x 2) 
vacancy structure,!^! our (\/3 x y/3) structure has only 
three-fold coordinated surface atoms and might thus be 
more favorable. 
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FIG. 4. Top and side views of the sample at 2350 K show- 
ing coexistence of the original Pandey chain (lower chains 
connected by solid lines) and the dimerized (alternating 
solid-dashed lines) reconstructions. The four atoms in the 
first bilayer connected by a dotted line in the top view and 
given as dark balls in the side view are lifted off (a precursor 
of surface graphitization) . The dark atoms underneath in the 
second bilayer have changed their coordination from four to 
three to form a new 7r-bonded chain. 

In conclusion, we have performed off-lattice GCMC 
simulations of the clean diamond (111) surface structure 
based on the Brenner potential. A spontaneous transi- 
tion from the ideal (1 x 1) to the stable Pandey (2 x 1) 
reconstruction is obtained. We also find metastable re- 
constructions very close in energy with strong dimeriza- 
tion of the lower atomic chain, which are shown to co- 
exist with the Pandey chain reconstruction at tempera- 
tures ~ 2000 K. Besides, we find a deep local minimum 
for the vacancy stabilized x \^3) structure. These 
meta-stable structures have a larger number of three-fold 
coordinated atoms at the surface. The absence of consen- 
sus on the structural details and electronic structure of 
the clean (111) surface might be related to these surface 
structures, which are peculiar of carbon and have never 
been considered so far. 

A.V.P. and A.F. would like to thank E. Vlieg, F. van 
Bouwelen, J.J. ter Meulen, W. van Enckevort, J. Scher- 
mer and B.I. Dunlap for useful discussions. D.P., F.E. 
and E.T. acknowledge support by MURST. 



* on leave from the Institute of Crystallography, 117333 
Moscow, Russia; e-mail petukhov@sci.kun.nl 

' present address: Max Planck Institut fur Festkoper- 
forschung, D-70569, Stuttgart, Germany. 

1 Y.G. Gogotsi, A. Kailer, K.G. Nickel, Nature, 401, 663 
(1999); A. De Vita et al, Nature 379, 523 (1996). 

2 F. Banhart, P.M. Ajayan, Nature, 382, 433 (1996). 

3 J. Tersoff, Phys. Rev. Lett. 61, 2879 (1988). 

4 E.C. Sowa, G.D. Kubiak, R.H. Stulen, J. Vac. Sci. Technol. 
A 6, 832 (1988). 

5 T.E. Derry, L. Smit and J.F. van der Veen, Surf. Sci. 167, 
502 (1986). 

6 W.J. Huisman et al, Surface Sci., 396, 241 (1998). 

7 W.J. Huisman, J.F. Peters, J.F. van der Veen, Surface Sci., 
396, 253 (1998). 

8 K.C. Pandey, Phys. Rev. Lett. 47, 1913 (1981). 

9 K.C. Pandey, Phys. Rev. B 25, 4338 (1982). 

A.V. Petukhov, A. Fasolino, D. Passerone, F. Ercolessi, 
Phys. Stat. Sol. (a) 174, 19 (1999). 

11 D.W. Brenner, Phys. Rev. B 42,9458 (1990). 

12 For the (001) surface see A.V. Petukhov, A. Fasolino, Phys. 
Stat. Sol. (a), to be published. 

13 F. Celestini, D. Passerone, F. Ercolessi and E. Tosatti, Surf. 
Sci. 402-404, 886 (1998). 

14 D. Vanderbilt, S.G. Louie, Phys. Rev. B 30, 6118 (1984). 

15 DR. Alfonso, D.A. Drabold, S.E. Ulloa, Phys. Rev. B 51, 
14669 (1995). 

16 G. Kern, J. Hafner, J. Furthmiiller, G. Kresse, Surface Sci., 
352-354, 745 (1996). 

17 A. Scholze, W.G. Schidt, F. Bechstedt, Phys. Rev. B 53, 
13725 (1996). 

18 R. Graupner et al, Phys. Rev. B 55, 10841 (1997). 

19 J.B. Cui, J. Ristein, L. Ley, Phys. Rev. B 59, 5847 (1999). 

20 S. Iarlori et al, Phys. Rev. Lett. 69, 2947 (1992). 

21 C. Kress, M. Fiedler, F. Bechstedt, Europhys. Lett. 28, 433 

(1994) . 

22 G. Jungnickel et al, Mat. Res. Soc. Symp. Proa, 383, 349 

(1995) ; Phys. Stat. Sol. (a) 154, 109 (1996). 
Th. Frauenheim et al, Phys. Rev. B 48, 18189 (1993). 
R. Graupner et al, Phys. Rev. B 57, 12397 (1998). 
D.J. Chadi, Phys. Rev. Lett. 26, 4762 (1982). 
R. Seiwatz, Surface Sci., 2, 473 (1964). 
M.P. Allen, D.J. Tildesley, Computer Simulation of Liquids 
(Clarendon Press, Oxford, 1987). 

R. Stumpf, P.M. Markus, Phys. Rev. B 47, R16016 (1993). 
F. Ancilotto et al, Phys. Rev. Lett. 65, 3148 (1990). 
A nimations, xyz files and many other det ails can be found 
at frittp : //www. sci .kun.nl/tvs/ carbon/. 



W.C. Fan et al, Phys. Rev. Lett., 62, 1516 (1989). 
F. Ancilotto, A. Selloni, E. Tosatti, Phys. Rev. B, 43, 
R14726 (1991). 

F. Bechstedt, W.G. Schmidt, A. Scholze, Europhys. Lett. 
35, 585 (1996). 



4 



